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With its central role in conventional BCS superconductivity, electron-phonon coupling has ap- 
peared to play a more subtle role in the phase diagram of the high temperature superconducting 
cuprates. The added complexity of the cuprates with potentially numerous competing phases includ- 
ing charge, spin, orbital, and lattice ordering, makes teasing out any unique phenomena challenging. 
In this review, we present our work using angle resolved photoemission spectroscopy (ARPES) to 
explore the role of the lattice and its effect on the valence band electronic structure in the cuprates. 
We provide an introduction to the ARPES technique and its unique ability to the probe the effect 
of bosonic renormalization (or "kink") on the near-Ei? band structure. Our survey begins with the 
establishment of the ubiquitous nodal cuprate kink leading to the way isotope substitution has shed 
a critical new perspective on the role and strength of electron-phonon coupling. We continue with 
recently published work on the connection between the phonon dispersion as seen with inelastic 
x-ray scattering (IXS) and the location of the kink as observed by ARPES near the nodal point. 
Finally, we present very recent and ongoing ARPES work examining how induced strain through 
chemical pressure provides a potentially promising avenue for understanding the broader role of the 
lattice to the superconducting phase and larger cuprate phase diagram. 

Keywords: photoemission, high temperature superconductivity, cuprates, phonon, electron-phonon coupling, 
migdal-eliashberg, inelastic x-ray scattering 
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I. FOREWORD 

The challenge of understanding the origin of high tem- 
perature superconductivity in the cuprates stems from 
the complicated interplay of differing orders and phe- 
nomena believed to exist. The goal of this article is 
to focus on one such phenomenon, the role of the lat- 
tice coupling to electronic states. Though historically 
significant in conventional superconductivity, it has only 
lately been receiving attention as a potentially impor- 
tant player in the physics of the cuprate phase diagram. 
Over the course of this article, we will be addressing the 
following areas: 1) How Angle Resolved Photoemission 
Spectroscopy (ARPES) can be used to probe and better 
understand electron self-energy effects. 2) A brief history 
of the ARPES "kink" seen in the cuprates and how both 
the energy scale it defines and its ubiquity in these sys- 



terns open the door to continued debate over low energy 
excitations. 3) How the cuprate isotope effect illuminates 
issues such as the role of phonons, the nature of the coher- 
ent and incoherent parts of the electronic dispersion, the 
limitations of current theory, and the subtle way differing 
competing orders may be interrelated within these sys- 
tems. 4) Recent work mapping phonon dispersion and re- 
lating it to ARPES data to underscore how phonon mode 
nesting may relate to the observed kink. 5) A survey of 
recent and ongoing work examining the role lattice strain 
may play in understanding electron-phonon coupling and 
potentially the larger phase diagram. Thus, our goal in 
this review is not to argue how much significance the 
lattice has to the development of high temperature su- 
perconductivity, but rather the ways in which we find it 
manifesting within these systems. 



II. SIGNATURE OF ELECTRON-PHONON 
COUPLING IN ARPES 

The use of ARPES to study electron-phonon coupling 
could be seen as the union of two different approaches 
to the study of phonons. First, and perhaps most obvi- 
ous, is the mapping of phonon dispersions using inelastic 
scattering techniques such as inelastic neutron or x-ray 
scattering. This momentum space perspective is gener- 
ally the most intuitive manner for understanding phonon 
modes within the lattice. Still, if we are seeking informa- 
tion about how electronic states interact with phonons, 
it is, at best, an indirect technique. Historically, tun- 
neling measurements such as those done on conventional 
superconductors, such as Pb [1 , have provided insight 
into how electron-phonon coupling directly affects elec- 
tronic states near the Fermi energy Ep. The previously 
unexpected features seen in the spectra were then able 
to be explained within a strong coupling form of Migdal- 
Eliashberg theory [2 . Still, to be able to have both the 
direct information of how phonon modes affect electronic 
states yet seen within a momentum space perspective, 
requires a different approach, an approach that ARPES 
is well-suited to offer. 



A. ARPES and A{k,uj) Analysis 

Over the last decade, ARPES has become a truly 
unique experimental probe with an ever growing number 
of publications in the field of correlated electronic sys- 
tems. Its central ability to directly probe the single parti- 
cle spectral function, A(k,C(;), makes its experimental in- 
sights highly sought after by condensed matter theorists. 
With angular resolution approaching 0.1°, a steadily im- 
proving energy resolution exceeding ImeV, as well as nu- 
merous experimental advances involving highly localized 
beam spots (Nano- ARPES), spin resolution, and laser- 
based pump-probe experiments, ARPES continues to be 
and will likely remain on the cutting edge of experimen- 



tal solid state physics. Because of the central role that it 
plays in the work described in this review, we will begin 
with a brief overview of the theory of ARPES, with an 
emphasis on the analytical techniques which are critical 
for studying systems such as the high Tc cuprates. 

It is customary to write the ARPES photocurrent in- 
tensity, /(k,^^), as 



/(k,cj) = M{k,uj)f{uj)A{k,uj) 



(1) 



where A(k, lo) is the crucial single particle spectral func- 
tion, i.e. the imaginary part of the single particle Green's 
function, G(k, uj) with k referring to the crystal momen- 
tum, while UJ is energy relative to the chemical potential. 
Modifying the spectral function, M(k,cj) is the matrix 
element associated with the transition from the initial to 
final electronic state which can be affected by such things 
as incident photon energy and polarization as well as the 
Brillouin zone (BZ) of the photoemitted electrons. Fi- 
nally, f{oj) is the Fermi-Dirac function indicating that 
only filled electronic states can be accessed. Because of 
the temperature scales used for data in this review, we 
will not distinguish between the chemical potential and 
the Fermi energy, Ep^ which should match at T = for 
conductors. Since ARPES measures the electron removal 
part of A(k,a;), we use high and low energy to refer to 
large and small negative uj value, respectively. (Addi- 
tionally, 'Binding energy'" and 'Energy' are often used 
for the same axis in figures, differing by a minus sign.) 
As a final point, one might find the contribution of the 
matrix element M(k, uj) in Eq. fl] a serious issue to an 
accurate interpretation of 74(k, cj) from /(k, cj). In prac- 
tice, the UJ dependence is small over an energy range of 
order 0.1 eV while the k dependence of M(k,cj), though 
important to consider, is reasonably understood by the 
ARPES cuprate community for the range of the data 
presented here. 

Of particular importance to our exploration of bosonic 
mode coupling is how electron self-energy effects appear 
in our ARPES analysis. This is nicely done by intro- 
ducing the electron proper self-energy i;(k,a;) = Rei;(k, 
uj) -\- i Iml](k, uj) which contains all the information on 
electron energy renormalization and lifetime. This leads 
to the Green's and spectral functions given in terms of 
the electron self energy i;(k, cj) 



G(k,cj) = 



1 



•e(k)-S(k,a;) 



(2) 



A{k,uj) = 



ImS(k,cj) 



7r{uj- e(k) - Rel](k, uj))'^ + ImS(k, uj)^ 



(3) 

where e(k) is the single electron band energy, often re- 
ferred to as the bare band structure. Finally, causality 
requires that ReS(k, cj) and Iml](k,cj) are connected to 
each other by the Kramers-Kronig relation. 

With advancements in the late 1990's, the unit infor- 
mation of an ARPES experiment consists of a two dimen- 
sional intensity map of binding energy and momentum 



along a "cut" though momentum space. These two di- 
mensional maps offer us two natural and complementary 
methods for analysis. First, one can hold the energy value 
of the electronic states studied fixed and observe the pho- 
toemission intensity as a function of momentum, a mo- 
mentum distribution curve (MDC). Similarly, one can fix 
the momentum space position and observe photoemission 
intensity as a function of energy at that momentum value, 
an energy distribution curve (EDC). These two methods 
constitute the core techniques for analysis of the spectral 
function ACk^uj) using ARPES. 

Within our review, "MDC analysis" refers to the 
method of fitting Lorentzian distributions to features 
in the MDCs as is commonly done in the field. This 
method of analysis has been very successful and can be 
understood based on some basic conditions, specifically 
the condition of "local" linearity in both the self energy 
S(k, cj) and e(k). In order for each MDC at a given en- 
ergy uj to be described with a Lorentzian function, both 
S(k, cj) and e(k) need to be linear within the narrow en- 
ergy and momentum range corresponding to the width 
of the analyzed peak. This condition is expected to gen- 
erally hold since both I](k,cc;) and e(k) are able to be 
expanded using simple Taylor expansions in the follow- 
ing way: 



S(k,a;) ^ ^{kp{uj),uj) + ^k{kp{uj),uj){k - kp{uj)) (4) 
e(k) ^ e{kp{uj)) + v{kp{uj)){k - kp{uj)) (5) 

where kp is the peak position of the MDC at uj and 
Y.k{kp{uj),uj) = [dJ:/dk]k=kpioj)' Plugging in these ex- 
pressions into Eq. [3| we obtain the following equations: 



A{k,uj) 



T{oj) 



(6) 



7r{k-kp{uj))^^r{uj)^ 
ReT^{kp{iju)^iju) = UJ — e{kp{uj)) (7) 

lm^{kp{uj),uj) =r{uj)[v{kp{uj)) ^J:k{kp{uj),uj)] (8) 

It is worth noting that we need not make the stan- 
dard assumption of momentum independence of i;(k, a;) 
for these results to be valid. This is important since mo- 
mentum dependence does exist for states away from the 
nodal cut (the diagonal direction to the Cu-0 bonds.) It 
is the last two equations which provide us with a precise 
meaning for kp{uj) and T{uj) as determined in the MDC 
analysis and, thus, our determination of S(k, cj) from 
ARPES. Eq. [7| demonstrates that a reasonable assump- 
tion for e(k) is needed to determine ReT^{kp{uj)^uj) while 
lm.Ti{kp{uj)^uj) presents the additional challenge of requir- 
ing the derivative Tik{kp{oj),uj). This is further compli- 
cated since though i;(k, cj) is a causal function for a fixed 
k value, T^{kp{uj) ^ u) is not. Thus, one cannot invoke 
the Kramers-Kronig relation to relate real and imaginary 
parts. Nevertheless, so long as these considerations are 
kept in mind to prevent over- interpret at ion, qualitatively 
r(co') and kp{uj) do offer access to the causal i;(k, u) since 
important structures such as the ARPES kink appear in 
both self-energies. 



Before the instrumental advances which pushed the 
unit information of ARPES towards a two dimensional 
map, one-dimensional data was taken, making EDC anal- 
ysis the more traditional method. Indeed, there are many 
advantages of this line of analysis: 1) Fixed momentum 
helps simplify the matrix element contribution to the 
photocurrent. 2) Momentum is a good quantum num- 
ber in a single crystal approximation making the EDC a 
more physical quantity, opening up spectral weight sum 
rules, as well as providing a clear physical meaning to 
the dispersion of EDC peaks. 3) In principle, an EDC 
analysis should be able to provide us with the causal 
I](k, cj) throughout in the entire two dimensional plane 
rather than a particular path determined by kp in the 
plane as with MDC analysis. However, EDC analysis 
is uniquely complicated by contributions from the Fermi 
function cutoff, /(cj), as well as both elastic and inelastic 
photoelectron background. This leads to a challenging 
lineshape to analyze in practice. Still, employing a sim- 
ilar method of Taylor expansion analysis as used earlier, 
we can expand the self-energy locally near the EDC peak 
yielding 



S(k,cj) ^ S(k,cjp(k)) + ^^{k,ujp{k)){uj-ujp{k)) (9) 

where i;^(k, cj) is the cj-partial derivative of i;(k,cj). 
Like before, we can insert these expressions into Eq. [3] 
getting the following relations in the vicinity of the peak. 



A(k,cj) 



Z(k) 



r(k,a;) 



(11) 



ReS(k,a;p(k)) = ujp{k) - e(k) 
ImS(k,a;p(k)) = r(k,a;)/Z(k) - 

lm^UKujp{k)){uj - ujp{k)) (12) 
Z(k) = l/(l-ReE^(k,c^^(k))) (13) 

When we compare these with our results for MDC 
analysis, the complementary nature of these two ap- 
proaches begins to appear. Unlike the Lorenzian line- 
shape of the MDCs, the EDC lineshape is modified by 
an asymmetry, which makes EDC analysis less favor- 
able for extracting self energy near Ep than MDC anal- 
ysis. However, the spectra at large uj is better ana- 
lyzed with EDCs thanks to the spectral sum-rule requir- 
ing ^(k,^^) -^ 1/cj, leading us to consequently expect 
Z(k) -^ 1 and S^(k, cj) -^ 0. This means that the por- 
tion of the spectral function which we associate with inco- 
herent excitations should begin approaching a Lorenzian 
lineshape and is better explored with EDCs, although 
inelastic background contributions at higher energy re- 
main important. In contrast, MDCs at higher energies 
begin to be affected by momentum dependent matrix ele- 
ment contributions as well as potential deviations of e(k) 
from a locally linear behavior. Thus, with both tools 
in our ARPES arsenal, we can undertake a more com- 
plete understanding of self-energy effects as they appear 
in A{k^uj). 
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FIG. 1. (a-c) ARPES data taken at 25 K on optimally doped 
Bi2212 superconductor (Tc = 92K), for a cut along the FY 
direction through the nodal point in momentum space (in- 
dicated in panel a inset), (d-f) The same sample and ori- 
entation, but taken nearer the Brillouin zone (BZ) edge (or 
antinodal point) in momentum space (indicated in panel d in- 
set). (a,d) Raw ARPES data taken with a color scale where 
intensity increases from pale yellow to green to blue to white 
to red. Here, blue (white) corresponds to 1/2 (3/4) of the 
maximum intensity. (b,e) Same data but as an "MDC map," 
where each MDC has been normalized so that its maximum 
and minimum intensities are 1 and 0, respectively. (c,f) Same 
data but now each EDC has been appropriately normalized 
to create an "EDC map." Thick black arrows indicate the en- 
ergy of the bosonic mode while the A is the superconducting 
gap. Energy resolution used here is ^-^ 15 meV. 



ducting gap is immediately apparent in the MDC "back- 
bending" observed in panel (e) near Ep within the gap 
energy scale. Evidence of this gap disappears for the 
nodal cut (panels a-c) as expected. 2) An abrupt de- 
viation in the electron dispersion around 70 meV below 
Ef (large black arrow) for both cuts. In both cuts, this 
corresponds to slower electron dispersion at lower energy 
while there is faster dispersion at higher energies above 
the 70meV energy scale. 3) Focusing particularly on the 
off- nodal cut, one sees evidence, even in the raw map, 
of an intensity decrease forming a local minimum at the 
70meV energy scale. This lineshape, further enhanced by 
the EDC map (panel f), is known as a "peak-dip- hump" 
and is associated with the presence of self-energy effects 
due to the coupling of electrons with a bosonic mode lead- 
ing to a redistribution of the spectral weight in the EDC 
spectra. 4) Although one would expect quasi-particles 
in a Landau Fermi liquid paradigm to become sharper 
(i.e. longer living) as one approaches Ei?, it is significant 
that the kink energy scale also marks a sudden change 
in coherence. One can make out from the panels, in par- 
ticular the EDC maps, abrupt changes in the MDC and 
EDC linewidths as one passes from states above and be- 
low the kink energy scale. To bring it all together, we 
can define the ARPES kink as an energy crossover sep- 
arating sharp, slowly dispersing, coherent states nearer 
Ejr from broader, quickly dispersing incoherent states at 
higher energy. 



The Nodal Kink 



B. Visualizing the Kink with ARPES 

Turning our attention to the physics of electron- 
phonon coupling in the superconducting cuprates, our 
prior discussion on how self-energy manifests in the 
ARPES spectral function points us towards the now well 
known "kink" feature. As Eqs.[7| and pT] quickly indicate, 
a sudden increase in the real part of S(k,co') at a partic- 
ular energy cj, would lead to a deviation of the measured 
peak from the single electron band structure e(k) at this 
energy scale. 

The result is seen in Fig. [l] which shows superconduct- 
ing phase data taken on the well-studied cuprate Bi2212 
at its optimal doping (Tc=92K). The two ARPES cuts 
are taken for states both at (panels a-c) and off (panels 
d-f) the nodal point. The different visualization meth- 
ods used for each cut are designed to enhance some key 
characteristics of the ARPES kink phenomenon prior to 
a more detailed, quantitative approach involving fittings. 
As labeled in the figure caption, the "MDC map" allows 
us to track the MDC dispersion and width. This is simi- 
larly true for EDCs in the "EDC map." The color scaling 
is chosen to give the peak maximum and half maximum 
distinct colors, red and blue respectively. 

From these maps, we can observe the following fea- 
tures: 1) The anisotropic d-wave nature of the supercon- 



With the initial discovery of the ARPES "kink" in the 
superconducting cuprate Bi2212[3], we have begun to de- 
velop a fuller picture of how low energy many-body effects 
manifest in the cuprates. On the heels of this discovery 
and the resulting debate, a systematic study regarding 
the origin of this kink discovered the feature's remarkable 
ubiquity across all cuprate families and dopings accessible 
by ARPES [4 . Fig. |2] summarizes these results particu- 
larly in double-layered Bi2212 and single-layered Bi2201 
and LSCO showing that the kink in the nodal direction 
exists all across these systems at essentially the same en- 
ergy, ~ 70 meV. One can give an estimate the coupling 
const ant. A, between the electrons and this bosonic mode 
by comparing the ratio of the group velocities above and 
below the kink energy, (1 + A) = VHigh/^Low [S] From 
this analysis, one finds evidence for a trend between dop- 
ing and the strength of the mode with an enhancement 
of S(k,co') as one tends to the underdoped side of the 
superconducting dome (panel f). 

Additionally significant is the continued existence of 
the kink below and above the superconducting transi- 
tion temperature (panel d-e), casting doubt on scenarios 
based on superconducting gap opening and particularly 
the magnetic mode scenario [6HTT]. Comparing the pho- 
toemission data with the neutron phonon energy at q 
= (tt, 0) (thick red arrow in panel a) and its dispersion 
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FIG. 2. (a-e) ARPES data at the nodal point showing the 
ubiquity of the ARPES dispersion kink as determined from 
MDC fittings over a variety of cuprates systems, dopings, and 
temperatures above and below Tc. The kink energy is indi- 
cated by the thick arrow and the momentum is rescaled so 
that k' is 1 at 170 meV binding energy, (f) Estimating the 
electron coupling constant A for the different samples as a 
function of their doping. Filled triangles, diamonds, squares, 
and circles are LSCO, Nd-LSCO, Bi2201, and Bi2212 in the 
first BZ, respectively. Open circles are Bi2212 in the second 
zone. Different shadings represent data from different exper- 
iments. Figures from Ref. [4^. 



(shaded area) [121 US! it was proposed [4 that the nodal 
kink results from coupling between quasi-particles and 
this zone boundary in-plane oxygen-stretching longitudi- 
nal optical (LO) phonon. Although this is the highest 
phonon mode contributing to the kink, quasi-particles 
are also coupled to other low energy phonon modes [14]. 
In favor of the electron-phonon coupled system is the 
drop of the quasi-particle width (Fig. 3| below the kink 
energy and the existence of a well-defined peak-dip- hump 
in the EDCs, a signature of an energy scale within the 
problem, persisting up to temperatures much higher than 
the superconducting critical temperature [15 . 
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FIG. 3. ARPES MDC dispersion data taken optimally doped 
Bi2Sr2CaCu208+<5 (Tc = 92 K) as discussed in Ref. [28]. 
Nodal point data (see inset) comparing (a) dispersion and (b- 
c) MDC full width half max (FWHM) showing little change 
in the energy of the ARPES kink with Tc. The MDC FWHM 
is related to the ImS(k,C(;). 
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FIG. 4. ARPES MDC dispersion data taken optimally doped 
Bi2Sr2CaCu208+<5 (Tc = 92 K) as discussed in Ref. Q^. Near 
antinodal point data (see inset) comparing (a) dispersion and 
(b-c) FWHM showing the ARPES kink energy is shifted to 
lower energy when passing above Tc. The MDC FWHM is 
related to the ImS(k, cj). 



D. The Near Antinodal Kink 

Many-body effects near the antinodal region of the BZ 
(Cu-0 bond direction) had been suspected for some time 
from earlier ARPES studies of the cuprates where ev- 
idence of the aforementioned peak-dip-hump lineshape 
was reported [16l |T7j . Although controversy has existed 
regarding the role of bilayer band splitting on the ob- 
served spectra in the double layered Bi2212 compounds, 
the presence of the peak-dip- hump lineshape was initially 
interpreted in terms of a magnetic phenomenon observed 
in YBCO and Bi2212 by inelastic neutron scattering [181 - 



|20] . With the resolution of the bilayer splitting [2TH23) , a 
low energy kink of approximately 40meV near the antin- 
odal region was reported for Bi2212 [24H26] . The dis- 
appearance of the kink above Tc and the decrease of its 
strength moving away from the antinodal region has led 
people to interpret the onset of this kink as coupling to 
collective magnetic excitations [24^-'26^, despite the ab- 
sence of these excitations for more heavily doped samples 

More recent studies [27', ^28^ have reported that the near 
antinode kink also persists above Tc, as seen in Fig. [4] 
However, the energy of this kink shifts towards higher 
energy, from 40meV to 70meV, upon entering the super- 



conducting state. This shift is consistent with the open- 
ing of a 30meV gap below Tc- The persistence of this 
energy scale above Tc can be clearly seen both in the 
dispersion (Fig. [4^) and in the MDC width (Fig. |4]3-c). 
This observation has led to a new interpretation of the 
antinodal kink in terms of electron-phonon coupling. In 
this case, it was proposed that the responsible phonon, 
with the right energy and momentum, is the Bi^ mode 

Still, it should be noted that spin fluctuations also ex- 
ist in the normal phase. Indeed they have been used 
within marginal Fermi liquid (MFL) theory to provide 
the anomalous self energy [29^ and have traditionally 
been called on to describe ARPES data in the normal 
phase just above T*. What is interesting is that although 
MFL theory can explain the ARPES dispersion in panels 
a of Figs. [3] and [4J it cannot explain the drop in linewidth 
seen in panels c of these figures, presenting difficulties to 
the original theory ^29^. Nevertheless, what we will find 
in the next section may suggest a profound connection 
between the physics of the lattice and spin. 

The results presented so far clearly suggest that the 
electron-lattice coupling could not be so easily neglected 
in any microscopic theory of cuprate superconductivity. 
As discussed in Section [TTbI the kink and its energy con- 
sistently indicate a sudden transition from sharp coherent 
electronic excitations into broader more incoherent ones. 
This makes understanding the origin of this phonon mode 
and its energy scale critical since it has such a substantial 
effect on low energy electronic states. 
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FIG. 5. (a-b) ARPES dispersions of a surface state for (a) H 
monolayer on W and (b) D monolayer on W. The dispersions 
are determined by peak positions from EDO fits and appear 
as circles while the lines serve as guides to the eye. (c) These 
two dispersions are compared where H = blue and D = red. 
Figure from |31j 



III. ARPES ISOTOPE EFFECT IN BI2212 

The role of the isotope effect (IE) in establishing the 
electron-phonon nature of Cooper pairing for traditional 
BCS superconductors is well known. However, when we 
consider the IE in the superconducting cuprates, its ef- 
fect on Tc is substantially less, leading researchers away 
from the electron-phonon paradigm of the BCS super- 
conductors. But with the ubiquitous cuprate kink seen 
by ARPES, the importance of the lattice returns to the 
forefront. Additionally from our discussion in |IIB[ the 
kink brings up questions about the relationship between 
the coherent peak (CP) seen at lower energies near Ep 
and the incoherent peak (IP) seen at higher energies. In- 
deed, hole doping affects the formation of these peaks 
differently [30^ with the CP strongly affected while the 
IP appears minimally changed. Should we be thinking of 
the CP and IP as different objects or fundamentally con- 
nected to each other? In this light, the kink energy scale, 
and thus the phonon mode responsible for it, becomes in- 
creasingly significant as the key crossover between these 
two domains. It is in light of such questions that we un- 
dertook our ARPES study of the IE to better understand 
the role of the lattice in these issues. 



A. Prelude - Isotope Effect in ARPES Dispersion 
of H/W 

Using APRES to explore isotope substituted samples 
presents us with an entirely new and wide open field of 
study. As of our work, the only other study in the lit- 
erature explores the surface state on W induced by H 
chemisorption [31]. Fig. Is] shows data taken for H on W 
(indicated in blue) compared to D on W (indicated in 
red). In spite of the broad peaks due to the instrumental 
resolution when compared to data on the cuprates, ex- 
tracting the EDC peak dispersion clearly shows two types 
of dispersions akin to that seen with the cuprate kink 
studies. Specifically, the slower low energy and faster 
high energy dispersions were understood in terms of the 
CP and IP, respectively. These results could be com- 
pared to predictions from the strong coupling form of 
Migdal-Eliashberg (ME) theory which, as discussed ear- 
lier, already explained the existence of peak-dip-hump 
features in tunneling spectroscopy as coupling to phonon 
modes. There was fair agreement with ME since the high 
energy linewidth (noted from panels a-b) as well as the 
kink energy position (dotted lines in panels a-b) were 
found to approximately scale as l/^/{M) where M is 
the mass of H or D. Additionally as expected from ME 
theory, the electron-phonon coupling A ~ 0.5, and the 
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FIG. 6. (a-b) MDC maps of the nodal point electronic states 
for cuts along the F-Y direction, (a) The ^^O sample and (b) 
the ^^O substituted sample with the horizontal arrows indi- 
cating the shift in ARPES kink energy with oxygen isotope. 
(c) The MDC dispersions determined from the ^^O, ^^O as 
well as a re-substituted ^^O samples for the cuts in (a-b). (d) 
Cartoon illustration of the kink shift in (c). (e) Real part 
of the electron self-energy, ReS(k,C(;), determined from the 
MDC dispersion using a linear approximation for the single 
electron bare band. As before, the ARPES kink position, de- 
fined by the peak in Rei;(k,c<;), is shifted to higher energy as 
indicated by the arrows. 



linewidth at high energy was ~ hcjp. Finally, panel c 
illustrates how the dispersions near the kink energy are 
affected by the isotope change, deviating most substan- 
tially near the kink energy while decreasing along both 
directions in energy, consistent with ME theory. 



B. Isotope Effect in ARPES Dispersion of Bi2212 

When we compare the ARPES IE seen in the 
H/W system to that measured on optimally doped 
Bi2Sr2CaCu208+(5, we find surprisingly different behav- 
ior. First, let's consider the behavior of the ARPES kink 
energy as summarized in Fig. [6J Panels a-b are MDC 
maps, as discussed in Fig.[l] from data taken at the nodal 
point in the F-Y direction. The two samples examined 
contain ^^O and ^^O in their Cu-0 planes which, as indi- 
cated by the horizontal arrow, already reveals a potential 
shift in the kink to lower energy with the substituted ^^O 
sample. 



Analyzing this more carefully, we turn to the MDC 
fitted dispersion for both isotopes plus a re-exchanged 
sample, ^^Oue-exch^ whereby we can take a studied ^^O 
sample and re-exchange ^^O back into the lattice. This 
provides us with a unique check on the IE. As indicated 
in panel c, we can observe, from the dispersions, a subtle 
shift in the kink energy between the isotopes of approxi- 
mately 5meV. We can additionally quantify the kink by 
estimating the bare single electron dispersion (using a 
linear approximation) and extracting the real part of the 
electron self-energy, Rel](k,cj), as described by Eq. 7\ 
The location of the peak in Rel](k, uj) corresponds to the 
kink energy and we similarly observe a shift in this peak 
with isotope change. Thus, the IE does have measure- 
able effect on the nodal ARPES kink energy, as later work 
would also confirm [32,, further establishing its phonon 
origin. 

A second aspect of panel c worth noting is the energy 
range where the IE is most pronounced. The maximum 
change in the dispersion occurs at higher energies, partic- 
ularly, beyond the kink energy, and is nearly non-existent 
at energies closer to Ep. This is particularly significant 
when we compare this result to the H/W work where the 
greatest deviation occurs near the kink energy. Fig. [5j 
We can explore this further by looking at the EDCs taken 
from this nodal cut, as presented in Fig.[7| Panel a shows 
the dispersion of the EDC peak from where it crosses Ep 
at kj7 as a sharp CP, to higher binding energy where it 
broadens and becomes the IP. Consistent with the MDC 
dispersions, we see very little change in the lineshape be- 
tween the two isotopes for the CP near Ejr. However, 
the IP at higher energy has a lineshape clearly affected 
by the IE. 

In light of this change at higher energy as compared to 
localization around the kink, a crucial question to ask is 
at what energy does the IE go away? Panel b of Fig. [7| 
attempts to address this issue by following the MDC dis- 
persion to even higher energies. Apparently, the IE seems 
to disappear around an energy scale of 2 to 3 times the 
antiferromagnetic coupling constant J (where J=4t^/U in 
the t-J model). This could suggest a profound intercon- 
nection between the effects of the lattice and spin on the 
electronic states in the superconducting cuprates. Fur- 
ther work is needed to better understand the connections 
between these phenomena. 

Up until now, we have focused entirely on the nodal 
point. Fig. |8] provides MDC dispersions for F-Y slices 
moving outward from the nodal point towards the antin- 
ode, both above and below Tg. There are a few important 
observations to make from this data. First, the kink en- 
ergy shows a subtle shift of approximately 5meV for all 
momentum cuts. Second, from comparing panel a with 
b, it appears that the magnitude of the IE may be, for 
all curves, diminished above Tc. Third, the IE remains 
relatively weak near the node while comparatively more 
pronounced near the antinode leading to a general cor- 
respondence between the kink strength. A, and the IE 
at high energy. Plotting this IE shift with respect to 
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FIG. 7. (a) EDCs taken on optimally doped 

Bi2Sr2CaCu2 08+5 samples with different oxygen isotopes at 
T = 25 K. The EDCs are from the same cut as in Fig. |6] 
and indicated by the inset. The sharp coherent peak (CP) 
near Ef and broader incoherent peak (IP) at higher energy 
are identified. The CP has nearly no isotope dependence 
while the IP has a more substantial IE, most strongly seen in 
curve 4. The small peak at Ep in curve 6 is the well-known 
superstructure (SS) replica of the main band. This figure 
is from Ref. 33 . (b) MDC dispersion taken to higher 
binding energy indicating the IE, more pronounced for the 
IP as in (a), disappears again above roughly 2-3 times the 
antiferromagnetic coupling energy, J. 
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FIG. 8. MDC dispersions from cuts parallel to F-Y taken off 
node towards the antinodal point, (a) Data taken in the su- 
perconducting phase (T = 25K). The inset shows the isotope 
energy shift vs the isotope- averaged superconducting gap, A. 
The isotope shift is measured at the momentum value where 
the isotope-averaged binding energy is 220 meV. The appar- 
ent linear correlation indicated by the dashed line is indepen- 
dent of the binding energy used, (b) MDC dispersions from 
the same cuts measured above Tc (T = 100 K). The inset 
illustrates the location of the cuts. Figures from Ref. [33 . 



the isotope averaged superconducting gap, gives a linear 
relationship seen in the inset of panel a. Finally, and 
potentially most surprising, there appears to be a sign 
change between the two dispersions as we transition from 
the node towards the antinode, which also appears both 
above and below Tc- 

This sign change is significant since one can examine 
its location in momentum space. When we plot these 
crossover points, surprisingly we find that they fall along 
a line defined in momentum space as qco =0.21 7r/a, 
where a is the lattice constant of the Cu02 plane. This 
is illustrated in Fig. [9] This particular wavevector (panel 
a) is in excellent agreement with the charge ordering 
wavevector seen in the far underdoped Bi2212 cuprate 
at low temperatures as explored by STM [34] (panel b), 
implying that the high energy part of the electronic struc- 
ture is strongly coupled to the order parameter, which is 
in turn strongly coupled to the lattice. 

To understand why at the qco line the isotope ef- 
fect changes sign, we used a simple charge density wave 
formation model, to show how an ordering mechanism 
can affect the quasiparticles dispersion at all energies. In 
panels c-d we present the opening of a gap in the dis- 
persion at the qco vector, due to a charge density wave 
formation. Based on the report that the pseudogap tem- 
perature is strongly isotope dependent and increases for 
the ^^O sample [35-37 , we assume that the magnitude 
of the gap is different between the two isotope samples. 



e.g. larger for the ^^O sample (panel d). This automati- 
cally leads to the appearance of the sign change at qco 
that migrates to lower energy as we move away from the 
node, as observed in the data. Although further studies 
are needed, we belive that the main reason why the pseu- 
dogap opening due to an ordering phenomena has never 
been observed in any ARPES experiment so far, is likely 
due to the short range nature of such ordering [38j |39j . 



This result has led us to consider a possible correla- 
tion between the charge or spin instability, with 4 ~ 5a 
periodicity, observed in the underdoped regime [34l [40] 
and the lattice effects relevant at optimal doping. In 
summary, these results suggest that the lattice degrees 
of freedom play a key role in the cuprates to "tip the 
balance" towards a particular ordered state. Simply put, 
the raw electron-phonon interaction may be small, but 
it can assist a certain kind of order through a cooper- 
ative enhancement of both the assisted order and the 
electron-lattice interaction |35l[36]. In particular, using 
a model which considers this electron-coupling boson as 
the critical source of charge order fluctuations [41], the 
sign change observed in our dispersion can be well repro- 
duced [39j . 
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FIG. 9. (a) Fermi surface of the upper half of the first BZ. 
Diagonal lines indicate cuts used in Fig. [S] while circles in- 
dicate the location in momentum space of the sign change 
crossover point for each of those cuts. These lie on a line 
indicated by the wavevector, qco and illustrated by the col- 
ored regions, (b) STM data independently determining this 
wavevector taken from Ref . [34 . (c) Cartoon illustrating how 
the charge ordering wavevector, qco, can open a gap at a 
binding energy where the electronic states are nestable, (d) 
Additional cartoon illustrating how the splitting, if slightly 
different in magnitude between the isotopes, can explain the 
observed sign change in the bands and its evolution as the 
dispersions intersect with qco- 



C. Isotope Effect in ARPES Width of Bi2212 

In addition to the shape of the electron dispersion, in- 
formation about linewidth, r(a;), can be extracted and 
provide critical information on the electronic states. As 
previously discussed, MDC analysis is not ideally suited 
to determining the linewidth at high energies, leaving us 
to examine the linewidth as it is obtained from EDC anal- 
ysis. Additionally, we have already observed a significant 
difference between states above and below the kink en- 



ergy, even prior to our discussion of IE. Thus, Fig. [TO] di- 
vides up the electronic states into those roughly between 
Ejr and the kink energy (0 to 70meV) and those beyond 
the kink energy (70meV to ~250meV). Taking the aver- 
age change in width of the EDC peaks shown in Fig. [7| 
between the isotopes for each of these two energy regions 
we find IE on the linewidth is very similar to its effect on 
the dispersion: 1) It is very small for the low energy co- 
herent electronic states, while much more significant for 
the higher energy incoherent states. 2) The magnitude of 
the IE is small at the node for these higher energy states, 
but it grows more substantial as one moves towards the 
antinode. 3) The effect is roughly linear with respect 
to the isotope averaged superconducting gap. A, as seen 
in Fig. |8] Yet, it does differ from the IE in the disper- 
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FIG. 10. Isotope effect on EDC widths determined from 
ARPES data taken at the slices indicated in the smaller panel 
(same as Fig.^ within the superconducting phase. The width 
change for each cut comes from an average over binding ener- 
gies below and above the kink energy: low (70meV to 0) and 
high (~250meV to 70meV) energy. These widths are plotted 
with respect to the isotope averaged superconducting gap. A, 
of each cut, showing little change at low energies but signifi- 
cant change above the kink energy. Figure is from Ref. [42]. 



sion since it lacks the sign change previously discussed. 
As will be addressed in depth later, the corresponding 
ME IE linewidth change is much smaller, about 2meV, 
making the trend in the high energy linewidth a serious 
failing of the theory for explaining the cuprate IE. 



D. Doping Dependence 

So far, the data shown has been on optimally doped 
Bi2212 with a hole doping determined by our ARPES 
experiment to be x = 0.16. But the question naturally 
arises, how is the IE on the ARPES data affected by a 
change in doping? Panel a of Fig. 11 shows angle inte- 



grated photoemission data obtained on two samples at 
optimal doping (x = 0.16) and one at slight over doping 
(x = 0.18). The effect is dramatic that given such a small 
change of only 2%, the IE, which normally shifts the IP 
position by ~ 30meV, is substantially reduced. This work 
was initially inspired by a separate study claiming not to 
see the IE in optimally doped samples of Bi2212 [43]. 
However, from superstructure analysis and studying the 
MDC dispersions, the samples discussed in Ref. [43 were 
actually slightly over doped. This doping dependence is 
intriguing but potentially puzzling given such sensitivity. 
Still, work in other correlated electronic systems, such as 
the manganites, demonstrates that such a small doping 
change does cause a qualitative change in the electronic 
structure [44 . Thus, the IE has a strong sensitivity to 
optimal doping and, therefore, there exists a rapid change 
in the electron-lattice interaction near optimal doping. 



E. Failure of Conventional Explanations 

Although the work of [43 had proposed that the IE 
was not apparent from ARPES on Bi2212, there are ad- 
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FIG. 11. (a) Angle integrated photoemission data from two 
sets of data at optimal doping (x = 0.16) and one set at a 
slight over-doping(x = 0.18). Data has been normalized to 
the area under curves for energies 0.5 and above, (b) MDCs 
taken near Ef along the F-Y nodal direction as seen from 
the inset. The peaks correspond to the intersection of the cut 
with the main Fermi surface band (M) as well as first-order 
superstructure reflections (S) and a second-order reflection 
(S'r)' Overlaid on top of the red curve is a thin grey line 
representing the ^^O curve shifted to model a 0.01 change in 
doping for comparison. Panel (a) from Ref. [45] . 



ditional explanations that one could invoke to explain 
the presence of this subtle effect which are important to 
address. As we will show, all of these conventional ex- 
planations turn out to be inadequate. Indeed, the strong 
temperature dependence by itself rules out all of the fol- 
lowing explanations as possible candidates. Still, this po- 
sition can be made stronger when we consider only the 
large low temperature IE in light of the following scenar- 
ios. 



1. Doping Issue 

As the previous section discussed, establishing the dop- 
ing of our samples, particularly that their optimal doping 
has remained unchanged after the isotope substitution, is 
critical to establishing the veracity of our claims. Doping 
level was preserved during the sample growth process by 
annealing the two samples (^^O and ^^O) in the exactly 
same environment save the oxygen gas. Yet, ARPES on 
the cuprates provides us with in-situ signatures of the 
doping level and its consistency. Most notably, one can 
use the Fermi surface size to quite precisely determine 
doping level, taking advantage of the well-known super- 
structure reflections of the main hole band structure. 

Using the nodal cut and making use particularly of 
the 2nd order reflection from the opposite side of the F 
point, panel b of Fig. [TT] shows the associated MDCs near 
Ef for three samples, including the re-substituted sam- 
ple. All three curves have good agreement with respect 
to peak positions. Of particular interest is the peak S[. 
which is due to the 2nd order superstructure Fermi sur- 
face replica from the opposite side of F. This means any 



doping change in the sample will affect the distance be- 
tween M and 6*^ twice as fast as the distance between 
M and F, while the distance between S and M remains 
fixed by the superstructure wavevector. This makes the 
distance between M and S^^ a sensitive measure of doping 
change. In fact, plotted on top of the red ^^O curve is 
a grey curve modeling a doping change of 0.01 based on 
a tight binding fit. From looking at how the ^^O peak 
positions compare to this, it is clear that the uncertainty 
in doping value is well below 0.01. 

However, one might be initially inclined to argue that 
the doping values between the two isotope samples are 
different based on the small difference in energy gap that 
has been previously reported [33] which would explain 
the observed IE. But this argument loses its plausibility 
for a few reasons. First, examining the two sets of data 
we report [33 , the difference in energy gap Aie — Aig 
actually differed in sign, with an ^ — 4 meV change for 
Figure 1 of Ref. [33^ and an ~ 5 meV change for Figures 
2,3 of Ref. 1^. Clearly, a empirically consistent change 
of gap is not obvious from our data. Secondly, this be- 
comes more evident over the many measurements (~ 20) 
we've done finding the Aie — Ai8 gap difference averaging 
to zero with a less than 1 meV difference, even while each 
individual value of Aie or Aig fluctuated by as much as 
5 meV. This was consistent with a typical uncertainty in 
the gap value from other sources at the time [46^. Thirdly, 
even if there were a consistently measured gap difference 
of 5 meV and this were taken to mean the doping val- 
ues were different, it still does not explain the observed 
crossover behavior in Fig. |8] More quantitatively, the 
doping change implied by a difference of 5 meV in gap 
{Ax = 0.017) is not sufficient to explain the large IE 
we've seen. With the associated doping change converted 
into shift in peak position at high energy, it corresponds 
to 5 meV for the nodal cut, and only 10 meV for the 
near-antinodal cut. So, these numbers not only have the 
same sign, but are off in magnitude by a factor of 3 to 
4. Thus, doping considerations do not offer a convincing 
explanation of the observed IE. 



2. Alignment issue 

Another conventional explanation for our observed IE 
could be sample alignment since even a small misalign- 
ment could create an apparent difference in the observed 
dispersion. A careful examination led us to focus on two 
particular issues: the exact location of the F point and 
the relative azimuthal orientation of samples with respect 
to each other. 

For the cuts 0-6 discussed here, we have found some 
evidence of deviations from our alignment based on the 
Fermi surfaces M, S and 5^ as described in panel b of 
These slightly different momentum paths are 



Fig. 11 



indicated in Fig. [12] to give a sense of magnitude. How- 
ever, we have found that, as panels b,c demonstrate, the 
expected difference in the dispersion, based again on our 
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FIG. 12. (a) Estimated azimuthal variation for cuts at and 
far off node indicated by the teal colored cuts relative to the 
red cut (F-Y). (b-c) Tight binding model for the nodal and off 
nodal cuts respectively representing the expected variation in 
dispersion due to the misalignment. 



tight binding fit, is none for cut and roughly about 4 
times too small for our cut 6 when compared to our mea- 
sured IE. It is also worth noting that the IE was repro- 
duced for data taken with analyzer slits rotated 45° with 
respect to this geometry, parallel to the MY direction [33]. 
This geometry has the advantage of being less sensitive 
to azimuthal tilting. Thus, alignment issues are not a 
likely explanation for the observed data. 



3. Static Lattice Issue 

A final concern that should be considered is that a 
static lattice effect may be responsible for the large IE 
since there are no high quality structural studies to rely 
on for insight. However, when we consider differences in 
crystal structure for isotope exchanged LSCO, they are 
only 0.1%. Given that a static structural effect is known 
to be more common in the LBCO and LSCO systems 
than the Bismuth based cuprates, one may reasonably 
assume that static lattice effects are significantly small in 
the Bi2212 cuprates and are an unlikely explanation for 
the observed IE. This is particularly strengthened since 
even if a static lattice effect were significant, it would 
need to be a particularly complicated static distortion 
which would make explaining the IE crossover difficult 
to accomplish. 




FIG. 13. (a-b) Comparison of ARPES EDCs for the three 
samples from panel c of Fig. [6] each normalized to same peak 
height. The cuts a and b for each panel respectively are indi- 
cated in the inset spanning the two panels along with the ap- 
proximate location along the cut in k-space for the EDCs. (c) 
Migdal-Eliashberg (ME) simulations for the expected change 
to (a) and (b) from ME for ^^O and ^^O, under-predicting 
the observed IE as described in text. Figures from Ref. [42] . 




-0.1 



0-0.2 -0.1 

k - kp (Ti/a) 







FIG. 14. (a-b) Comparison of ARPES MDC dispersions for 
these two locations with the different isotopes, (c) Migdal- 
Eliashberg (ME) simulation for change in dispersion with iso- 
tope again under-predicting the magnitude and location of 
the expected band structure change. 



and theoretical [561458] works indicating an interaction 
strength beyond this theory. So, in light of our work, we 
will distinguish between those aspects which clearly go 
beyond ME theory as well as those more consistent with 
the theory. 

To illustrate this former point, Fig.fTs) begins by show- 



F. Beyond the Migdal-Eliashberg Picture 

With the IE results more soundly established in light 
of other potential explanations, we return our attention 
to the ME theory and ask whether its applicability is still 
appropriate for the results seen. One would initially ex- 
pect ME theory to offer the best theoretical model of the 
observed ARPES kink in this doping region of the phase 
diagram. However, in view of the discussion at the start 
of this section, it is not entirely obvious that ME the- 
ory can be used to describe the broader incoherent spec- 
tral weight. This is a significant concern because of the 
prevalent use of ME theory in the context of the ARPES 
kink in the cuprates despite many experimental 
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FIG. 15. (a) Simulation of EDCs at k = ki? for the off nodal 
cut seen in Fig. [l3j) modeled with a small polaron theory 
[59] . (b) The expected broadening of this lineshape showing 
the small IE at low energy (see panel a inset) with a more 
substantial IE shift at higher energy. Figures from Ref. [42] . 
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ing EDCs for two cuts, a nodal(panel a) and a near antin- 
odal (panel b). Although subtle near the nodal point, 
there is a significant deviation between the ^^O (or the 
re-substituted ^^Oi^) and the ""^^O which is not effectively 
modeled within the expected change from ME theory 
(panel c) and already mentioned previously in regards 
to the IE on ARPES width. This certainly comes as no 
surprise since we already knew the behavior of the dis- 
persion near the kink energy for the cuprate data was not 
well modeled by ME as compared to the H/W data which 
better follows its predictions. This is again emphasized 
in Fig. [l4j where the IE on the MDC dispersions near 
the node and antinode (panels a-b) is much bigger when 
compared to the substantially smaller expected change 
from the ME theory (panel c) as well as the aforemen- 
tioned localization of the change to just near the kink 



energy. Additionally, as again illustrated in Fig. [Mj the 
presence of the momentum dependent sign change in the 
dispersion (panels a-b) would defy explanation by any 
simple application of ME theory at least in the absence 
of an additional ordering phenomenon. 

Fundamentally, the failure of the ME theory has its 
origin in the single phonon loop approximation for the 
electron self-energy, resulting in a small IE, particularly 
at higher energy. This suggests to us that we need to in- 
crease the electron-phonon coupling in our model. To ac- 
complish this, we employ a Holstein model in the strong- 
coupling limit as described in Ref. [59]. The resulting 
comparison of this small polaron theory to the data can 
be found in Fig. 15 These simulate the lineshape at k 



= ki? for the off nodal cut b. The multiple peaks seen 
in panel a occur due to the strong multi-phonon "shake- 
up" peaks which appear at harmonics of ujp. With the 
expected broadening due to the phonon continuum and 
strong electron-electron interactions, the result, modeled 
in panel b, successfully reproduces the weaker IE for the 
low energy CP as well as the larger IE for the broader 
IP (^30). Additionally, it produces an ARPES linewidth 
which is more realistic. 

These results clearly indicate that the ME theory is 
insufficient for describing the ARPES data. However, it 
should be noted that this strong coupling theory used is 
not quite the right solution either. Within the strong cou- 
pling theory, the IP is not expected to have any dispersion 
while the CP is expected to be nearly non-dispersive as 
well. Furthermore, it predicts a very small quasi-particle 
weight (Z<< 0.1) for the CP (panels c-d) which is not ob- 
served. Both of these issues are better modeled by ME 
theory, so it is important to note that these shortcom- 
ings in the theory are overcome by weakening the inter- 
action strength. This leads us to propose that the proper 
paradigm for understanding self-energy in the optimally 
doped cuprates is an intermediate regime. In this regime, 
there is a significant multi-phonon contribution to elec- 
tron self-energy and both the CP and IP show strong dis- 
persions. Though one may have reasonably expected the 
IE to also weaken with diminished coupling, some studies 
[6Q| r61j show anomalous enhancement of the IE at inter- 



mediate couplings. Nevertheless, the key message here is 
that as far as modeling the electron-phonon physics seen 
in IE, the ME theory is not sufficient and that any subse- 
quent theory must incorporate electron-phonon coupling 
with a strength beyond the ME paradigm. 



IV. LINKING IXS AND ARPES THROUGH 
PHONONS 

With this renewed interest regarding the role of the lat- 
tice in the superconducting cuprates as seen by ARPES, 
a natural direction to explore is to more carefully map 
the phonon dispersion in addition to that of the electron. 
A comparison of these dispersions can reveal particularly 
interesting physics such as when a phonon wavevector 
matches 2ki?, where kp is the Fermi momentum, lead- 
ing to the well-known Kohn anomaly. Additional, nest- 
ing of the Fermi surface in systems with particularly 
strong electron-lattice coupling can drive the formation 
of charge density waves. Thus, it becomes important to 
compare data that directly probes the phonon mode dis- 
persions within the Cu-0 plane, such as one can get from 
inelastic neutron scattering (INS) and x-ray scattering 
(IXS) as discussed in Section [ill Then we can compare 
these results with the observations of the ARPES kink 
particularly in light of the potential change in the binding 
energy of the kink as one moves from the nodal point (60- 
70meV) towards the antinodal point (30-40meV). Fur- 
ther establishing the phonon nature of the nodal kink as 
well as shedding light on the lower energy antinodal kink 
is something that a combined ARPES and IXS study 
could achieve. 



A. IXS Measurements on La-Bi2201 

We turn our attention to the single layered 
Bi2Sri.6La0.4Cu2O6+5 (La-Bi2201) cuprate which has 
several advantages for such a study. Like other Bismuth 
superconducting cuprates, the sample surfaces are high 
quality for ARPES experiments. These samples have 
never shown any evidence of magnetic resonance modes, 
simplifying the comparison between ARPES kink and 
scattering by removing a potentially additional bosonic 
mode to couple with the electrons. Moreover, no experi- 
mental reports of the optical phonon dispersion exist on 
these materials to date. The challenge for scattering is 
the lack of large single crystals, effectively ruling out an 
INS experiment. Additionally, though IXS can probe 
the sub- millimeter single crystals available, there is a 
very low inelastic cross section associated with the bond 
stretching (BS) mode, a likely candidate for the nodal 
kink, making observing the mode challenging. Still, both 
IXS and INS have observed evidence in the past of the 
Cu-0 bond stretching (BS) mode at the metal-insulator 
phase transition in the superconducting cuprates. 

Fig. [16] encapsulates the IXS experiment with panels a 
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FIG. 16. (a-b) Raw IXS phonon spectrum taken from op- 
timally doped La-Bi2201 at lOK. (a) Focusing on the high 
energy part and the LO phonon dispersions for Q= (3+^,0,0) 
with ^ ranging from the BZ center to the BZ face, (b) Stronger 
low energy part of the phonon spectrum. Solid lines indicate 
fits, dashed lines show the elastic tail, and dotted lines in- 
dicate the modes, (c-d) The peaks of these dispersions are 
plotted with cosine dashed lines for the two potential dis- 
persion scenarios: crossing (c) and anti-crossing (d) Figures 
taken from [62] . 
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FIG. 18. (a) MDC dispersions measured for three different 
momentum cuts along the F-Y orientation with cut 1 at the 
nodal point while cuts 2 and 3 are further toward the BZ 
boundary, near to the edge of the pseudogap phase Fermi 
arcs [65]. (b) Experimentally determined Fermi surface with 
the cuts from (a) indicated. The solid line indicates a con- 
stant energy contour at the kink energy, 63meV, while the 
shadow area indicates the region where the nodal kink ap- 
pears bounded by the indicated nesting wavevectors. The 
inset shows the IXS dispersion and peak FWHM (seen as er- 
ror bars) of the BS mode. Note: The apparent Fermi arcs 
seen are due to the experimental resolution. Figures can be 
found in |62J. 



zone center but they disperse, becoming indistinguish- 
able around ^ of 0.22-0.25. When ^=0.45, the two peaks 
clearly emerge again, leading to two potential scenarios 
of panels c and d depending on the symmetry of the two 
branches. If they have the same symmetry, they anti- 
cross (panel d), otherwise they simply cross (panel c). 
Our attempts at distinguishing between these two sce- 
narios using classical shell model calculations could not 
reproduce the low and high energy modes observed in a 
reliable way. Thus, we have been unable to distinguish 
between the potential scenarios. Finally, we compared 
our data with other IXS data in the literature, summa- 
rized in Fig. [T7| We find broad agreement for a softening 
of the BS mode (panel a) between ^=0.2-0.3 as well as a 
maxima in the full-width half-max (FWHM) of the BS 
mode peak. 



FIG. 17. (a) Bond stretching mode softening in a variety 
of optimally hole doped cuprates compared to the La-Bi2201 
data, (b) Similar comparison of peak FWHM. Figures from 
[62j . Data on LSCO and LaBaCuO come from ^63^ and the 
Hgl201 comes from [64] 



and b illustrating the relative weakness of the BS phonon 
peak relative to both the elastic line as well as other 
modes. Focusing on these higher energy longitudinal op- 
tical modes, we map out their dispersion across the BZ 
from the center (^ = 0) towards the BZ face. Panel 
a shows this evolution where the red peak is identified 
as the BS mode and the results of which are plotted in 
panels c and d. We see the two distinct peaks at the 



B. ARPES Measurements on La-Bi2201 

These results take on an deeper meaning when we com- 
pare them with our ARPES studies on La-Bi2201. Fig. [18 
displays the ARPES results in comparison to the IXS 
data. Panel a shows MDC analysis of the electronic dis- 
persions taken at the nodal point (curve 1) and away from 
the node (curves 2 and 3) as indicated by the slices along 
the Fermi surface in panel b. We see the well established 
higher energy kink at 63±5meV for the nodal cut. As 
we move away from the node, this kink abruptly disap- 
pears between curves 2 and 3, replaced with only a lower 
energy kink of 35meV. It is significant that this shift oc- 
curs at the tips of the so-called "Fermi Arc," region of the 
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Fermi surface which becomes ungapped at temperatures 
above the superconducting Tc but below the so-called 
pseudogap temperature, T* [66]. Beyond the Fermi arc 
as one moves towards the BZ edge, the gap reopens again 
for reasons which remain mysterious. It is this transition 
point in momentum space between the arc and where the 
gap opens that curves 2 and 3 straddle, though our data 
was taken in the superconducting phase. 

When we compare this with the IXS data as seen in the 
inset of panel b, we discover that the 63meV kink has an 
energy that corresponds well to the softened BS mode. 
Even more interesting, as the grey shaded region in panel 
b illustrates, the region where the 63meV kink is observed 
corresponds to a section of the Fermi surface nestable by 
wavevectors within the softened part of the phonon mode 
dispersion from IXS. The sudden disappearance of this 
kink between curves 2 and 3 corresponds to the stiffening 
of the BS mode when ^ < 0.22. A final insight is that 
this BS mode is supposed to be non-dispersive at about 
80meV along other momentum directions, in particular 
the [110]. We see no strong feature above 63meV in our 
data, helping confirm that the nodal charge carriers are 
preferentially coupled to the softer Cu-0 half-breathing 
BS mode which disperses along the [100] direction, a re- 
sult suggested by local spin-density approximation + U 
results [67]. This work not only provides additional di- 
rect evidence for the lattice origin for the 60-70meV kink 
seen near the node, it associates it with the softened Cu- 
O half-breathing BS mode along the [100] direction. Its 
ability to nest the Fermi surface topology once again un- 
derscores the importance of electron-phonon interactions 
to the physics of the superconducting cuprates. 
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FIG. 19. (a) Superconducting Tc for optimally doped 
Bi2Sri.6Ln0.4Cu2O6 for a series of substituted Lan- 
thanides (Ln) with increasing atomic radius mismatch, 
AR= I i?sr - RLn\. See Ref. [72]. (b) Cartoon illustrating the 
location of the substituted Lanthanide right above the Cu-0 
plane, (c) Nodal point EDCs illustrating the quasi-particle 
peak for samples with increasing strain, colored in (a), from 
La to Eu. Inset quantifies the half-width half-max of the 
peaks for these samples. 



[TT] . Thus, better understanding of this aspect of the role 
of the lattice is important to our general understanding 
of electron-lattice physics in these systems. 



V. LATTICE STRAIN IN BI2201 

Our most recent work has continued this exploration 
into the role of the lattice from yet another perspective. 
There has been growing independent work from a vari- 
ety of experimental probes J68H7QJ positioning the role of 
the lattice not simply as a source of self-energy effects 
on the near-Ei? low energy electronic states but poten- 
tially as an additional axis within the hole-doped phase 
diagram. Specifically, it is the effect of lattice strain, 
both external and internal via chemical pressure, which 
offers us this new axis to the cuprate phase diagram af- 
fecting the superconducting dome. Work using external 
pressure has indicated critical pressures where the Tc ap- 
pears to saturate for a range of cuprate hole dopings [68^ 
as well as being coupled to other physical quantities sug- 
gesting a significant new critical point along this axis 
[69^. With chemical pressure, work has suggested that 
combining doping with strain on the Cu-0 layer also re- 
veals that the true quantum critical point is shifted along 
the strain axis [70]. Additionally, effects related to lat- 
tice disorder, particularly in the Sr-0 blocks nearest the 
Cu-0 planes, may also have a dramatic effect on the for- 
mation of the superconducting phase within the cuprates 



A. Lanthanide Substituted Bi2201 

Due to experimental considerations, the best method 
for introducing strain into the lattice for an ARPES study 
is via chemical pressure. Specifically, we can use Lan- 
thanide substituted single-layered Bi2Sri.6Lno.4Cu06 to 
access this strain in a tunable way [72]. All the sam- 
ples were grown at optimal doping to simplify the anal- 
ysis, making the focus solely on the tuning parameter of 



strain. As Fig. 19 illustrates, the substitution of the Lan- 
thanide elements for the Strontium right above the crit- 
ical Cu-0 plane (panel b) leads to a monotonic decrease 
in the measured Tc of the samples. The essential variable 
to quantify this Tc-competing strain is the atomic radius 
mismatch, AR, as seen on the abscissa of panel a, which is 
determined by the difference between the Strontium and 
the substituted Lanthanide atomic radii, \Rsr — ^Ln|- 

We have been able to take data on samples throughout 
this spectrum of radius mismatch. Although these sam- 
ples allow us a lattice-based, tunable parameter which 
competes with superconductivity, one can pose the ques- 
tion if this should be thought of within a lattice strain 
or lattice disorder paradigm. Panel c of Fig. 19 pro- 
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FIG. 20. (a) Nodal point MDC dispersions taken from T-Y 
cuts on four different samples of increasing strain (La, Pr, 
Nd, and Eu). Lines serve as guides to determine deviation 
from the expected dispersion. Horizontal shaded regions cor- 
respond to the two potential ARPES kink energy scales, (b- 
c) Estimating the electron coupling A on the electronic states 
from each of the two regions indicated in (a), (bl-2) Higher 
energy kink A for each strain as a function of (1) Lattice mis- 
match A and (2) Tc. (cl-2) Lower energy kink A analyzed 
like (b). 



vides evidence that, at least for the nodal states, the 
strain paradigm appears valid. With increasing lattice 
mismatch, the width of the quasi-particle peak does not 
increase but, on close examination, may even be decreas- 
ing with increased strain. The introduction of lattice dis- 
order should decrease the lifetime of the electronic states, 
seen as an increase in the width of the CP. It also has 
been broadly suggested within the ARPES cuprate com- 
munity that observing a sharp quasi-particle peak at the 
nodal point is necessary for confirming that the cleaved 
surface can provide trustworthy ARPES results. 



B. ARPES Kink vs. Strain 

Throughout this review, we have continually discussed 
the ARPES kink, particularly near the node, in terms of 
electron-phonon coupling. So considering the expected 
effect of the substituted lanthanides on lattice strain and 
the view that these states at least can be understood 
within a strain paradigm, we focus our attention again on 



the MDC dispersions at the nodal point. Fig. 20 presents 
our findings for substituted Ln=La, Pr, Nd, Eu [73] . 
There are two aspects of these dispersions we wish to 



focus on in the context of this review. First, in agreement 
with the earlier work on La-Bi2201 discussed in Section 
IV Bl we can observe a kink around 55-60meV which re- 



mains at that energy, for the most part, throughout the 
strain spectrum. However, what appears to change with 
strain is the electron-coupling constant, A associated with 
the renormalization of these states. In the same man- 
ner as was done for the cuprate systems described in 
Section [II C[ we can estimate A for this mode, which is 
plotted in panel b. We find the strength of this mode 
appears enhanced by the increasing strain of the lattice 
mismatch with a generally linear behavior. Equivalently, 
one can plot A as a function of sample Tc (panel b inset) 
and, as one would expect, there is a negative, linear re- 
lationship between the superconducting Tc and strength 
of this phonon mode. Secondly, we find that although 
it appears linear for the La-Bi2201 at energies less than 
the 60meV kink, the more strained compounds appear 
to have additional rounding of the band structure nearer 
to Ei?, most obvious in the highly strained Eu-Bi2201. 
This mode appears to be around 25-30meV which could 
be important since this is closer to the mode energy ob- 
served near the antinodal point (as discussed in Section 
II D) from the peak-dip-hump EDC lineshape. We have 



also observed this from the kink in the MDC analysis 
beyond the nestable region of the IXS softened phonon 
mode in Section llVBI 

As with the higher energy mode, we can attempt to 
extract the A from this more elusive mode independent 
of the 60meV kink, the result of which is seen in panel 
c. Unlike the higher energy mode which is apparent in 
all strain, this lower energy feature appears to turn on 
at the node only as strain is introduced, leading to a 
broadly linear relationship similar to the higher energy 
mode. The origin of this feature, remains mysterious, po- 
tentially related to an apical oxygen mode particularly in 
light of the location of the substituted lanthanide above 
the Cu-0 plane. 

These results are significant for at least three reasons. 
First, it has been a prevailing thought that the electronic 
states at the nodal point are uniquely unaware of the 
entry into the superconducting phase. With the d-wave 
symmetry of the gap function, the nodal point states are 
the only electronic states which still cross through Ep 
with no gap opening. Additionally, one could argue that 
the continued appearance of a sharp quasi-particle at the 
node is merely because these states are protected from 
the superconducting physics. However, clearly the affect 
of lattice strain can be seen on these states and the weak- 
ening of the superconducting state is present in the elec- 
tronic dispersion of the nodal quasi-particles. Secondly, 
one finds still additional evidence that the ~60meV kink 
has its origin in the physics of the lattice. Even more 
intriguing, the appearance and potential enhancement of 
the lower energy kink with lattice strain would tenta- 
tively suggest its origin also is somehow connected to the 
lattice (such as the aforementioned apical oxygen mode) 
and not merely a magnetic mode at the nodal point. It 
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could be pointed out that the lanthanides do carry with 
them magnetic moments whose effect on the Cu-0 plane 
is far from understood. While La has no magnetic mo- 
ment, Pr, Nd, and Eu all have experimentally determined 
magnetic moments of around 3.5/iB [74 • Whether this 
may be related to the sudden appearance of this mode at 
the node for samples beyond La requires further study. 
Finally, at least for the 60meV mode, one finds evidence 
that this phonon mode is somehow connected to the for- 
mation of the superconducting phase. From its behavior, 
it appears to be related to a competing order, associated 
with the lattice, which may be affecting the formation of 
the superfiuid. 



VI. SUMMARY 



to the ARPES kink near the nodal point. From this, we 
find the electronic states within the BZ can be divided up 
accordingly, with those nearest to the node (often argued 
to be responsible for the superfiuid) uniquely impacted by 
the mode wavevector's ability to nest them. Finally, the 
potential for understanding the role of the lattice via the 
additional axis of pressure/strain has already produced 
remarkable results for the evolution of the nodal elec- 
tronic states and their self-energy along this axis as well 
as pointing towards additional insights. In all, our work 
has and continues to explore how the cuprate electronic 
structure is affected by electron-phonon coupling as yet 
another pillar in our nearly 25 year quest to construct a 
full understanding of this unconventional superconduct- 
ing phase. 



Throughout this work, the re-occurring theme has 
been the growing importance that the lattice and its cou- 
pling with electronic states has within the still mysteri- 
ous phase diagram of the hole doped cuprates. Beginning 
with the ubiquitous nodal kink and its likely origin in the 
bond stretching phonon mode, we have explored the ef- 
fect of the lattice on both the coherent and incoherent 
parts of the near Ep electronic band structure through 
the IE. This has confirmed that, near optimal doping, 
the traditional ME model must give way to a stronger 
electron-phonon coupling, yet still shy of a polaronic pic- 
ture. We have traced out the dispersion of the BS phonon 
mode with IXS and have found its important connection 
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